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2016  FINAL  REPORT  -  JEFFREY  A.  MAGEE,  M.D.,  Ph.D. 

1.  INTRODUCTION 

The  goal  of  this  project  is  to  understand  why  FLT3-lnternal  Tandem  Duplication  ( FLT3-ITD )  mutations 
cause  acute  myeloid  leukemia  (AML)  more  frequently  in  adults  than  in  children.  FLT3-ITD  encodes  a 
constitutively  active  FLT3  tyrosine  kinase.  This  mutation  occurs  in  ~30%  of  adult  AML,  but  only  5- 
10%  of  pediatric  AML  and  <1%  of  infant  AML.  This  suggests  that  the  mutation  preferentially 
transforms  older  progenitors.  To  test  whether  FLT3-ITD  has  age-specific  effects  on  hematopoietic 
stem  cells  (HSCs)  and  other  hematopoietic  progenitors,  we  characterized  its  function  in  fetal, 
neonatal  and  adult  mouse  FISCs.  We  discovered  that  FLT3-ITD  selectively  alters  the  fate  of  adult, 
but  not  fetal  progenitors,  primarily  due  to  normal  temporal  changes  in  transcriptional  regulation  that 
make  adult  progenitors  competent  to  express  FLT3-ITD  target  genes. 

2.  KEYWORDS 

•  Flt3-lnternal  Tandem  Duplication  (FLT3-ITD) 

•  Hematopoietic  stem  cell  (HSC) 

•  Acute  myeloid  leukemia 

•  Stat5 

•  MAP-kinase 

3.  ACCOMPLISHMENTS 
Major  goals: 

Aim  1:  To  test  whether  FLT3-ITD  depletes  HSCs,  expands  restricted  progenitors  and  promotes  a 
myeloproliferative  neoplasm  during  the  adult,  but  not  fetal  stage  of  development. 

Aim  2:  To  test  whether  fetal  and  adult  hematopoietic  progenitors  have  different  FLT3-ITD  driven 
signal  transduction  mechanisms  and  gene  expression. 

Aim  3:  To  test  whether  ectopic  Lin28b  expression  impedes  FLT3-ITD  driven  HSC  depletion  and 
leukemogenesis. 

Accomplishments  under  these  goals: 

Aim  1:  To  test  whether  FLT3-ITD  depletes  HSCs,  expands  restricted  progenitors  and  promotes  a 
myeloproliferative  neoplasm  during  the  adult,  but  not  fetal  stage  of  development. 

We  completed  all  of  the  proposed  experiments,  and  the  data  are  incorporated  into  a  manuscript  that 
was  just  accepted  for  publication  at  eLife  (acceptance  date  11/21/2016).  Please  see  attached 
manuscript. 


Aim  2:  To  test  whether  fetal  and  adult  hematopoietic  progenitors  have  different  FLT3-ITD  driven 
signal  transduction  mechanisms  and  gene  expression. 

We  completed  all  of  the  proposed  experiments,  and  the  data  are  incorporated  into  a  manuscript  that 
was  just  accepted  for  publication  at  eLife  (acceptance  date  11/21/2016).  Please  see  attached 
manuscript. 
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Aim  3:  To  test  whether  ectopic  Lin28b  expression  impedes  FLT3-ITD  driven  HSC  depletion  and 
leukemogenesis. 

We  have  completed  all  of  the  proposed  experiments.  Ectopic  Lin28b  expression  did  not  prevent  Flt3- 
ITD  driven  HSC  depletion  or  leukemogenesis.  We  are  now  evaluating  other  fetal  gene  products,  and 
we  are  dissecting  the  gene  regulatory  networks  that  make  adult  hematopoietic  progenitors  sensitive 
to  FLT3-ITD. 

Opportunities  for  training  and  professional  development: 

This  work  was  selected  for  an  oral  presentation  at  the  American  Society  of  Hematology  annual 
meeting  (December  2015),  and  it  has  formed  the  basis  for  an  R01  proposal  that  just  received  a  5th 
percentile  score  (final  council  review  pending).  Thus,  the  work  supported  by  CDMRP  has  opened  new 
lines  of  investigation,  and  funding  opportunities,  that  will  sustain  my  lab  for  years  to  come.  I  received 
guidance  from  my  primary  mentor,  Sean  Morrison,  throughout  the  funding  period. 

Dissemination  of  results  to  communities  of  interest: 

This  work  has  been  presented  in  oral  presentations  at  the  International  Society  for  Stem  Cell 
Research  Annual  meeting  and  the  American  Society  of  Hematology  annual  meeting.  As  noted  above, 
an  manuscript  describing  our  results  has  just  been  accepted  at  eLife. 

Plans  for  the  next  reporting  period  to  accomplish  the  goals: 

Nothing  to  report. 

4.  IMPACT 

Impact  on  the  development  of  the  principal  discipline  of  the  project: 

This  project  provides  the  first  evidence  that  leukemia  causing  mutations  can  have  age-specific  effects 
on  blood  forming  stem  cells  and  other  immature  blood  cells.  This  likely  contributes  to  the  differences 
between  pediatric  and  adult  leukemias,  and  It  raises  the  possibility  of  targeting  developmental 
programs  to  suppress  leukemogenesis. 

Impact  on  other  disciplines: 

Nothing  to  report. 

Impact  on  technology  transfer: 

Nothing  to  report. 

Impact  on  society  beyond  science  and  technology: 

Nothing  to  report. 
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5.  CHANGES/PROBLEMS: 


Changes  in  approach/reasons: 

Nothing  to  report. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them: 

Nothing  to  report. 

Changes  that  had  a  significant  impact  on  expenditures: 

Nothing  to  report. 

Significant  changes  in  the  use  or  care  of  human  subjects: 

Nothing  to  report. 

Significant  changes  in  the  use  or  care  of  vertebrate  animals: 

Nothing  to  report. 

Significant  changes  in  the  use  of  biohazards  and/or  select  agents: 

Nothing  to  report. 

6.  PRODUCTS 
Publications: 

Portern  SN,  Cluster  AS,  Yang  W,  Busken  KA,  Patel  RM,  Ryoo  J,  Magee  JA.  Fetal  and  neonatal 
hematopoietic  progenitors  are  functionally  and  transcriptionally  resistant  to  Flt3-ITD  mutations.  eLife 
[in  press]. 

Books  or  one-time  publications: 

Nothing  to  report. 

Other  publications: 

Nothing  to  report. 
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Websites: 


Nothing  to  report. 


Technologies  and  techniques: 

Nothing  to  report. 

Inventions,  patent  applications  or  licenses: 

Nothing  to  report. 

Other  products: 

Nothing  to  report. 


7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
Individuals  who  have  worked  on  this  project: 


Name: 

Jeffrey  Magee,  M.D.,  Ph.D. 

Project  role 

PI 

Nearest 
person  month 
worked 

7.2 

Contribution 
to  Project 

Dr.  Magee  oversaw  the  entire  research  effort  including  the 
design,  conduct  and  interpretation  of  experiments.  He  also 
performs  experiments. 

Funding 
Support 
besides  DOD 

St.  Baldrick's  Foundation,  Hyundai  Hope  on  Wheels, 
Gabrielle’s  Angel  Foundation 

Name: 

Shaina  Porter,  Ph.D. 

Project  role 

Post-doctoral  fellow 

Nearest 
person  month 
worked 

12 

Contribution 
to  Project 

Dr.  Porter  designed,  conducted  and  interprets  experiments 
under  the  guidance  of  the  PI.  She  has  conducted  the 
experiments  related  to  temporal  changes  in  FLT3-ITD 
regulation  of  HSCs  and  the  signal  transduction  and  gene 
expression  studies. 

Funding 

Support 

Hyundai  Hope  on  Wheels 
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Name: 

Andrew  Cluster,  M.D. 
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Nearest 
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to  Project 

Dr.  Cluster  characterized  the  FLT3-ITD;  Runxl  compound 
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Funding 
Support 
besides  DOD 

Departmental  support  for  clinical  trainees  -  graduate  medical 
education  support. 

Name: 

Jenna  Voigtmann 

Project  role 

Research  technician 
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worked 

6 
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to  Project 

Ms.  Voigtmann  managed  the  mouse  colony  and  assisted  with 
experiments. 

Funding 
Support 
besides  DOD 

St.  Baldrick’s  Foundation  and  institutional  funds. 
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Kelsey  Busken 

Project  role 

Research  technician 

Nearest 
person  month 
worked 
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Support 
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Name: 
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Research  technician 
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Ms.  Patel  co-managed  the  mouse  colony  and  assisted  with 
experiments. 

Funding 
Support 
besides  DOD 
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Name: 

Wei  Yang 

Project  role 

Bioinformaticist 

Nearest 
person  month 
worked 

0.5 

Contribution 
to  Project 

Dr.  Wang  performed  all  of  the  computational  work  associated 
with  this  project  -  e.g.  microarray  analyses. 

Funding 
Support 
besides  DOD 

No  grant  support  -  supported  by  the  department  of 
biostatistics 

Name: 

Jiyeon  Ryoo 

Project  role 

Research  technician 

Nearest 
person  month 
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4 

Contribution 
to  Project 
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experiments. 

Funding 
Support 
besides  DOD 

Hyundai  Hope  on  Wheels,  St.  Baldrick’s  Foundation 

Changes  in  active  other  support  of  the  PD/PUI  or  senior  key  personnel: 

Nothing  to  report. 


Other  organizations  that  were  involved  as  partners: 

Nothing  to  report. 


8.  SPECIAL  REPORTING  REQUIREMENTS 

Nothing  to  report. 


9.  APPENDICES 

Please  see  attached  manuscript 
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29  Abstract 

30  The  FLT3  Internal  Tandem  Duplication  ( FLT3itd )  mutation  is  common  in 

31  adult  acute  myeloid  leukemia  (AML)  but  rare  in  early  childhood  AML.  It  is  not 

32  clear  why  this  difference  occurs.  Here  we  show  that  Flt3ITD  and  cooperating 

33  Flt3ITD/Runx1  mutations  cause  hematopoietic  stem  cell  depletion  and  myeloid 

34  progenitor  expansion  during  adult  but  not  fetal  stages  of  murine  development.  In 

35  adult  progenitors,  FLT3itd  simultaneously  induces  self-renewal  and  myeloid 

36  commitment  programs  via  STAT5-dependent  and  STAT5-independent 

37  mechanisms,  respectively.  While  FLT3itd  can  activate  STAT5  signal  transduction 

38  prior  to  birth,  this  signaling  does  not  alter  gene  expression  until  hematopoietic 

39  progenitors  transition  from  fetal  to  adult  transcriptional  states.  Cooperative 

40  interactions  between  Flt3ITD  and  Runxl  mutations  are  also  blunted  in 

41  fetal/neonatal  progenitors.  Fetal/neonatal  progenitors  may  therefore  be  protected 

42  from  leukemic  transformation  because  they  are  not  competent  to  express 

43  FLT3itd  target  genes.  Changes  in  the  transcriptional  states  of  developing 

44  hematopoietic  progenitors  may  generally  shape  the  mutation  spectra  of  human 

45  leukemias. 

46 
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47  Introduction 

48  Acute  myeloid  leukemia  (AML)  can  occur  at  any  stage  of  life  yet  the 

49  mutations  that  cause  AML  differ  between  childhood  and  adulthood,  especially 

50  when  one  compares  young  children  to  adults  (Chaudhury  et  al.,  2015).  For 

51  example,  MLL  translocations  and  GATA1  mutations  are  common  in  infant  and 

52  early  childhood  AML  but  rare  in  adult  AML  (Andersson  et  al.,  2015;  Horton  et  al., 

53  2013;  Pine  et  al.,  2007).  Mutations  in  FLT3,  NPM1,  DNMT3A,  TET2  and  IDH1 

54  are  all  common  in  adult  AML  but  rare  in  infant  and  early  childhood  AML  (Cancer 

55  Genome  Atlas  Research  Network,  2013;  Ho  et  al.,  2011;  Liang  et  al.,  2013; 

56  Zwaan  et  al.,  2003).  The  genetic  differences  between  pediatric  and  adult  AML 

57  are  not  absolute,  but  they  reflect  a  more  general  phenomenon  in  leukemia 

58  biology  -  leukemias  in  infants,  young  children,  older  children  and  adults  have 

59  different  genetic  and  epigenetic  landscapes,  different  mechanisms  of 

60  transformation  and  different  clinical  courses  (Downing  and  Shannon,  2002). 

61  Efforts  to  interpret  AML  genomes  and  translate  the  information  into  useful 

62  therapies  will  need  to  account  for  the  influences  of  age  and  developmental 

63  context  on  leukemia  cell  biology.  This  will  require  a  better  understanding  of  how 

64  normal  developmental  programs  shape  the  process  of  leukemogenesis. 

65  The  mutations  that  cause  AML  are  thought  to  accrue  first  in  pre-leukemic 

66  hematopoietic  stem  cells  (HSCs)  or  committed  hematopoietic  progenitor  cells 

67  (HPCs)  (Jan  et  al.,  2012;  Welch  et  al.,  2012),  and  several  properties  of  these 

68  cells  change  between  fetal  and  adult  stages  of  life:  1 )  Fetal  HSCs  divide 

69  frequently  and  retain  their  self-renewal  capacity  through  cumulative  division 
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70  cycles  (Pietras  and  Passegue,  2013).  In  contrast,  adult  HSCs  are  usually 

71  quiescent,  and  self-renewal  capacity  declines  with  cumulative  divisions  (Foudi  et 

72  al.,  2009;  Pietras  and  Passegue,  2013;  Wilson  et  al.,  2008).  2)  Fetal  and  adult 

73  HSCs  have  distinct  self-renewal  mechanisms.  For  example,  Soxl  7  is  required  for 

74  fetal,  but  not  adult,  HSC  self-renewal  (Kim  et  al.,  2007).  Etv6,  Ashll,  Mil  and  Pten 

75  are  all  required  for  adult,  but  not  fetal,  HSC  self-renewal  (Hock  et  al.,  2004; 

76  Jones  et  al.,  2015;  Jude  et  al.,  2007;  Magee  et  al.,  2012).  3)  Fetal  and  adult 

77  HSCs  give  rise  to  committed  progenitors  with  distinct  epigenetic  landscapes 

78  (Huang  et  al.,  2016;  Xu  et  al.,  2012)  and  distinct  lineage  biases  (Benz  et  al., 

79  2012;  Copley  et  al.,  2013;  Yuan  et  al.,  2012).  These  observations  raise  the 

80  question  of  whether  mutations  can  have  age-specific  effects  on  gene  expression, 

81  self-renewal,  differentiation  and  ultimately  leukemogenesis.  If  so,  competence  for 

82  transformation  may  be  a  heterochronic  property  of  HSCs  and  HPCs,  and  this 

83  may  explain  why  pediatric  and  adult  leukemias  have  different  mutations. 

84  The  FLT3  Internal  Tandem  Duplication  ( FLT3itd )  is  an  example  of  an  AML 

85  driver  mutation  that  occurs  more  commonly  in  adults  than  in  young  children  (30- 

86  40%  of  adult  AML,  5-10%  of  AML  in  children  <10  years  old,  <1%  of  infant  AML) 

87  (Meshinchi  et  al.,  2006).  FLT3itd encodes  a  constitutively  active  tyrosine  kinase 

88  receptor  that  has  been  shown  to  activate  the  STAT5,  MAP-kinase  (MAPK),  PI3- 

89  kinase  (PI3K),  STAT3  and  NF-kB  signal  transduction  pathways  in  various 

90  contexts  (Choudhary  et  al.,  2007;  Gerloff  et  al.,  2015;  Radomska  et  al.,  2006). 

91  Mice  with  a  targeted  Flt3ITD  mutation  develop  myeloproliferative  neoplasms 

92  (MPN)  (Lee  et  al.,  2007;  Li  et  al.,  2008),  and  several  other  mutations  (e.g.  Npml, 
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93  Tet2  and  Runxl  mutations)  cooperate  with  Flt3ITD  to  drive  AML  in  mice  much  as 

94  in  humans  (Mead  et  al.,  2013;  Mupo  et  al.,  2013;  Rau  et  al.,  2014;  Shih  et  al., 

95  2015).  In  the  absence  of  cooperating  mutations,  Flt3ITD  drives  adult  HSCs  into 

96  cycle  and  depletes  the  HSC  pool  (Chu  et  al.,  2012).  This  may  explain  why 

97  FLT3itd  mutations  occur  late  in  the  clonal  evolution  of  human  AML  -  adult  HSCs 

98  must  first  acquire  mutations  that  preserve  (or  ectopically  establish)  self-renewal 

99  capacity  in  pre-leukemic  progenitors  -  but  it  also  raises  the  question  of  why 

100  fetal/neonatal  HSCs,  which  have  an  inherently  high  self-renewal  capacity  (He  et 

101  al.,  2009),  do  not  give  rise  to  FLT3itd  positive  AML  more  often  than  is  observed. 

102  To  better  understand  how  developmental  context  shapes  myeloid 

103  leukemogenesis,  we  characterized  the  effects  of  Flt3ITD  on  HSC  self-renewal, 

104  myelopoiesis,  signal  transduction  and  gene  expression  at  several  stages  of  pre- 

105  and  post-natal  development.  Flt3ITD  did  not  cause  HSC  depletion  or  myeloid 

106  progenitor  expansion  until  after  birth.  This  was  true  even  in  the  presence  of  a 

107  cooperating  Runxl  loss-of-fu notion  mutation.  The  FLT3itd  protein  phosphorylated 

108  STAT5  during  both  pre-  and  post-natal  stages  of  development  while  it  hyper- 

109  activated  the  MAPK  pathway  only  after  birth.  To  our  surprise,  MAPK  inhibition 

110  failed  to  rescue  HSC  depletion  and  myeloid  progenitor  expansion  in  adult  Flt3ITD 

111  mice,  and  Stat5a/b  deletion  greatly  exacerbated  these  phenotypes.  FLT3itd 

112  target  genes,  including  STAT5  targets,  were  not  induced  in  fetal  HSCs  or  HPCs 

113  despite  pre-natal  STAT5  phosphorylation.  Instead,  FLT3itd  target  gene  activation 

114  coincided  with  a  normal  transition  from  fetal  to  adult  gene  expression  that  was 
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115  evident  by  two  weeks  after  birth.  These  temporal  changes  in  FLT3itd  target  gene 

116  expression  were  observed  even  in  the  setting  of  a  cooperating  Runxl  mutation. 

117  Our  data  establish  a  crucial  role  for  developmental  context  in  the 

118  pathogenesis  of  FLT3itd- driven  AML.  Fetal  and  neonatal  progenitors  are 

119  protected  from  transformation  because  they  are  not  competent  to  express 

120  FLT3itd  target  genes.  This  likely  explains  why  FLT3itd  mutations  are  more 

121  common  in  adults  than  young  children,  and  it  may  reflect  a  more  general  role  for 

122  developmental  programming  in  leukemia  pathogenesis. 

123 

124  Results 

125  Flt3ITD  does  not  deplete  fetal  FISCs 

126  Since  FLT3itd  occurs  more  commonly  in  adult  AML  patients  than  in  young 

127  children,  we  hypothesized  that  it  might  have  age-specific  effects  on  self-renewal 

128  and  myelopoiesis.  We  first  tested  whether  Flt3  expression  changes  with  age.  We 

129  measured  Flt3  transcript  expression  in  CD150+CD48'Lineage'Sca1+c-kit+  HSCs 

130  and  CD48+Lineage'Sca1+c-kit+  HPCs  from  8-10  week  old  adult  and  embryonic 

131  day  (E)14.5  fetal  mice  by  quantitative  RT-PCR  (qRT-PCR).  Flt3  was  more  highly 

132  expressed  in  HPCs  than  in  HSCs  at  both  ages  (Fig.  1  A),  consistent  with  prior 

133  studies  (Buza-Vidas  et  al.,  201 1 ),  but  its  expression  did  not  change  with  age  in 

134  either  cell  population  (Fig.  1A).  Flow  cytometry  confirmed  that  the  FLT3  protein  is 

135  expressed  in  both  fetal  and  adult  progenitors  (Fig.  IB). 

136  Since  Flt3ITD  has  previously  been  shown  to  deplete  adult  HSCs  (Chu  et 

137  al.,  2012),  we  tested  whether  the  mutation  has  a  similar  effect  on  fetal  HSC 
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138  numbers.  We  measured  HSC  numbers  in  8-10  week  old  adult  bone  marrow  and 

139  El 4.5  fetal  livers  from  wild  type,  Flt3ITD/+  and  Flt3ITD/ITD  mice.  Adult  Flt3ITD/+  mice 

140  had  ~50%  fewer  HSCs  than  wild  type  littermates,  consistent  with  prior  studies, 

141  and  Flt3ITD/ITD  mice  had  a  near-complete  loss  of  phenotypic  HSCs  (Fig.  1C).  HSC 

142  depletion  in  the  bone  marrow  was  not  accompanied  by  extramedullary  expansion 

143  of  HSCs  in  the  spleen  (Fig.  1 E),  in  contrast  to  other  leukemogenic  mutations  (e.g. 

144  Pten  deletion)  that  cause  depletion  of  bone  marrow  HSCs  but  marked  expansion 

145  of  the  spleen  HSC  population  (Magee  et  al.,  2012;  Porter  et  al.,  2016).  Unlike 

146  adult  mice,  Flt3ITD/+  and  Flt3ITD/ITD fetal  mice  had  similar  numbers  of  HSCs  as  wild 

147  type  littermates  (Fig.  1 D).  The  Flt3ITD  mutation  therefore  depletes  adult,  but  not 

148  fetal  HSCs. 

149  We  next  tested  whether  fetal  Flt3ITD/ITD  HSCs  are  functionally  impaired.  We 

150  performed  limiting  dilution  transplantation  assays  with  either  8-10  week  old  adult 

151  bone  marrow  cells  (600,000,  100,000,  50,000  or  10,000  CD45.2  donor  cells 

152  competed  with  300,000  CD45.1  adult  bone  marrow  cells)  or  El 4.5  fetal  liver  cells 

153  (100,000,  50,000  or  10,000  CD45.2  donor  cells  competed  with  300,000  CD45.1 

154  adult  bone  marrow  cells).  Two  independent  experiments  were  performed,  and 

155  fetal  and  adult  donor  cells  were  transplanted  at  the  same  time  in  each 

156  experiment.  Multi-lineage  reconstitution  was  assessed  every  4  weeks  for  16 

157  weeks  following  the  transplants,  and  functional  HSC  frequencies  were  calculated 

158  by  Extreme  Limiting  Dilution  Analysis  (Hu  and  Smyth,  2009).  Adult  Flt3ITD/ITD  bone 

159  marrow  had  significantly  fewer  functional  HSCs  than  adult  wild  type  bone  marrow 
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160  (P<0.00001,  Fig.  IF).  In  contrast,  wild  type  and  Flt3ITD/ITD  fetal  livers  had  similar 

161  HSC  frequencies  (Fig.  1G). 

162  Our  findings  raised  the  question  of  whether  fetal  Flt3ITD/ITD  HSCs  can 

163  mature  and  become  depleted  after  transplantation  into  adult  recipient  mice.  To 

164  test  this,  we  measured  donor  HSC  chimerism  in  primary  recipients  of  100,000 

165  wild  type  and  Flt3ITD/ITD  fetal  liver  cells  (from  Figure  1 G).  Donor  Flt3ITD/ITD  HSCs 

166  were  significantly  depleted  in  primary  recipient  mice,  but  wild  type  competitor 

167  HSCs  were  not  (Fig.  1 H).  Overall  donor  bone  marrow  chimerism  was  not 

168  significantly  different  between  recipients  of  wild  type  and  Flt3ITD/ITD  fetal  liver  cells 

169  (Fig.  1 1).  Secondary  transplants  confirmed  depletion  of  Flt3ITD/ITD  HSCs  in  the 

170  marrow  of  primary  recipients  (Fig.  1 J,  K).  Thus,  fetal  Flt3ITD/ITD  HSCs  are 

171  functional,  but  they  lose  repopulating  activity  after  transplantation  into  adult 

172  recipient  mice. 

173 

174  Flt3ITD  causes  FISC  depletion  and  myeloid  progenitor  expansion  by  2  weeks  after 

175  birth  in  mice 

176  We  next  sought  to  define  the  age  at  which  Flt3ITD  begins  to  deplete  HSCs 

177  and  expand  myeloid  progenitor  populations.  We  measured  HSCs,  HPCs  and 

178  granulocyte-monocyte  progenitor  (GMP)  numbers  at  E14.5,  El 6.5,  post-natal 

179  day  (P)0  and  P14.  Flt3ITD/+  and  Flt3ITD/ITD  mice  had  normal  HSC  numbers  at  all 

180  ages  prior  to  birth  (Fig.  2A).  HSC  depletion  was  evident  at  P14  in  both  the  bone 

181  marrow  and  the  spleen,  though  not  to  the  extent  observed  in  adult  bone  marrow 

182  (Fig.  2B).  We  observed  a  modest  increase  in  Flt3ITD/+  and  Flt3ITD/ITD  HPCs  and 
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183  GMPs  at  PO  (Fig.  2C,  E),  and  this  phenotype  became  more  severe,  particularly  in 

184  Flt3ITD/ITD  mice,  by  P14  (Fig.  2D,  E).  Spleen  enlargement  due  to  MPN  was  evident 

185  by  P14  in  Flt3ITD/+  and  Flt3ITD/ITD  mice,  but  E14.5,  El  6.5  and  PO  liver  sizes  were 

186  not  increased  relative  to  wild  type  littermates  (Fig.  2F).  These  data  show  that 

187  Flt3ITD  causes  HSC  depletion,  HPC/GMP  expansion  and  MPN  beginning  at  or 

188  shortly  after  birth. 

189  FLT3itd  mutations  usually  occur  late  during  the  clonal  evolution  of  human 

190  AML.  This  raises  the  question  of  whether  fetal/neonatal  Flt3ITD  mice  can  exhibit 

191  HSC  depletion  and  HPC/GMP  expansion  when  a  cooperating  mutation  is 

192  present.  To  test  this,  we  analyzed  HSC,  HPC  and  GMP  frequencies  in  Flt3ITD/+; 

193  Runx1f/+;  Vavl-Cre  mice.  Both  mono-  and  bi-allelic  RUNX1  loss-of-fu notion 

194  mutations  co-occur  with  FLT3itd  in  human  AML  (Schnittger  et  al.,  201 1 ),  and 

195  Runxl  deletions  synergize  with  Flt3ITD  to  cause  AML  in  mice  (Mead  et  al.,  2013). 

196  For  the  purposes  of  these  studies  we  focused  on  mono-allelic  Runxl  deletions 

197  because  bi-allelic  deletions  severely  depleted  phenotypic  HSCs  irrespective  of 

198  the  Flt3  genotype  (data  not  shown).  These  effects  were  likely  due  to  previously 

199  described,  Runxl- dependent  changes  in  CD48  expression  (Cai  et  al.,  2011). 

200  We  evaluated  HSC,  HPC  and  GMP  frequencies  in  1)  Flt3+/+;  Runx1f/f  or 

201  Runx1f/+;  Cre-negative  (control),  2)  Flt3ITD/+;  Runx1f/f  or  Runx1f/+;  Cre-negative 

202  ( Flt3ITD/+ ),  3)  Flt3+/+;Runx1f/+;  Vavl-Cre  ( Runx1A/+ )  and  4)  Flt3ITD/+;Runx1f/+;  Vavl- 

203  Cre  ( Flt3ITD/+;  Runx1A/+)  littermates  at  E14.5,  P0,  P14  and  P21 .  HSCs  were 

204  severely  depleted  in  P14  and  P21  Flt3ITD/+:  Runx1^+  mice  relative  to  controls  and 

205  single  mutant  mice  (Fig.  3C,  D).  In  contrast,  all  four  genotypes  of  mice  had 
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206  similar  HSC  frequencies  at  P0  (Fig.  3B),  and  Runxl  heterozygosity  increased 

207  HSC  frequency  at  El 4.5  irrespective  of  the  Flt3  genotype  (Fig.  3A).  HPCs  and 

208  GMPs  were  markedly  expanded  in  P14  and  P21  Flt3ITD/+:  Runxl^  mice  (Fig.  3G, 

209  H,  K,  L).  These  populations  were  only  modestly  expanded  in  compound  mutant 

210  mice  at  P0  (Fig.  3F,  J),  and  they  were  not  expanded  at  all  at  E14.5  (Fig.  3E,  I). 

211  We  next  tested  whether  compound  Flt3ITD  and  Runxl  mutations  had  age- 

212  specific  effects  on  HSC/HPC  function.  We  transplanted  100,000  P0  liver  cells  or 

213  P21  bone  marrow  cells  from  control  or  Flt3ITD/+:  Runx1N+  littermate  donors,  along 

214  with  300,000  wild  type  competitor  cells,  into  irradiated  CD45.1  recipient  mice.  At 

215  2  weeks  after  the  transplants,  we  observed  CD45.2+  donor-derived  leukocytes  in 

216  the  peripheral  blood  of  all  recipients,  irrespective  of  donor  age  or  genotype  (Fig. 

217  3M,  N).  At  4  weeks  after  the  transplants,  donor  chimerism  was  significantly  and 

218  dramatically  reduced  in  recipients  of  Flt3ITD/+;  Runx1A/+  P21  donor  cells  as 

219  compared  to  recipients  of  control  P21  donor  cells  and  Flt3ITD/+:  Runx1A/+  P0  donor 

220  cells  (Fig.  3M).  Indeed,  only  1  of  15  recipients  of  Flt3ITD/+:  Runx1A/+  P21  donor 

221  cells  had  multi-lineage  donor  chimerism  (>0.5%  CD45.2+  myeloid  and  lymphoid 

222  cells)  (Fig.  30).  In  contrast,  all  recipients  of  control  and  Flt3ITD/+:  Runx1A/+  P0 

223  donor  cells  had  multi-lineage  donor  chimerism  (Fig.  30).  These  differences  were 

224  evident  even  when  we  focused  specifically  on  myeloid  chimerism  (Fig.  3N,  P),  so 

225  they  were  not  simply  a  reflection  of  altered  lineage  biases  in  the  Flt3ITD/+: 

226  Runx1N+  progenitors.  Altogether,  these  data  show  that  Flt3ITD  has  developmental 

227  context-specific  effects  on  HSC  depletion,  myeloid  progenitor  expansion  and 

228  repopulating  activity,  even  when  paired  with  a  cooperating  Runxl  mutation. 
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229 

230  Flt3ITD  activates  STAT5  and  MAPK  signal  transduction  pathways  in  adult  HSCs 

231  and  HPCs,  yet  it  only  activates  STA  T5  in  fetal  HSCs 

232  To  better  understand  why  Flt3ITD  has  developmental  context-specific 

233  effects  on  HSCs  and  HPCs,  we  sought  to  better  characterize  the  pathways  that 

234  mediate  FLT3itd  signal  transduction  in  vivo.  We  isolated  25,000  HSC/multipotent 

235  progenitors  (HSC/MPPs;  CD48'Lineage'Sca1+c-kit+),  HPCs  and  GMPs  from  adult 

236  mice  by  flow  cytometry,  and  we  performed  Western  blots  to  assess 

237  phosphorylation  of  STAT5,  STAT3,  ERK1/2  (a  MAPK  pathway  protein)  and  AKT 

238  (a  PI3K  pathway  protein).  Both  STAT5  and  ERK1/2  were  hyper-phosphorylated 

239  in  Flt3ITD  HSC/MPPs  and  HPCs,  as  well  as  in  adult  Flt3ITD;  Runxl**  AML  cells 

240  (Fig.  4A-C).  In  contrast,  STAT3  and  AKT  were  not  hyper-phosphorylated  in 

241  Flt3ITD  mutant  HSC/MPPs  or  HPCs  (Fig.  4A  and  Fig.  4  -  figure  supplement  1  A), 

242  and  Rictor  deletion  (PI3K/mTORC2  pathway  inactivation)  did  not  rescue  Flt3ITD- 

243  driven  HSC  depletion  or  MPN  (Fig.  4  -  figure  supplement  1 B,  C).  These  findings 

244  suggest  that  the  STAT5  and  MAPK  pathways  mediate  FLT3itd  signal 

245  transduction  in  hematopoietic  progenitors,  but  the  STAT3  and  PI3K  pathways  do 

246  not. 

247  We  next  tested  whether  FLT3itd  signal  transduction  changes  between 

248  fetal  and  adult  stages  of  development.  We  isolated  HSC/MPPs  and  HPCs  from 

249  wild  type  and  Flt3ITD  mice  at  E14.5,  P0,  P14  and  8  weeks  after  birth.  We 

250  performed  Western  blots  to  assess  STAT5  and  ERK1/2  phosphorylation.  STAT5 

251  was  hyper-phosphorylated  in  Flt3ITD  mutant  HSC/MPPs  and  HPCs  at  all  stages 
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252  of  development,  though  the  degree  of  STAT5  phosphorylation  appeared  to 

253  increase  with  age  (Fig.  4D,  E).  ERK1/2  was  only  hyper-phosphorylated  in  post- 

254  natal  Flt3ITD  mutant  HSC/MPPs  and  HPCs  (Fig.  4D,  F).  Several  other  signal 

255  transduction  proteins,  including  STAT3,  AKT,  ribosomal  protein  S6,  p38  and 

256  JNK,  were  not  hyper-phosphorylated  in  Flt3ITD  HSC/MPPs  or  HPCs  at  any  age 

257  tested,  or  their  phosphorylation  was  undetectable  (data  not  shown).  Our  data 

258  reinforce  other  studies  that  have  implicated  STAT5  and  MAPK  as  key 

259  downstream  effectors  of  FLT3itd  signaling  (Choudhary  et  al.,  2007;  Radomska  et 

260  al.,  2006).  However,  the  data  suggest  that  these  pathways  are  not  coupled  - 

261  STAT5  is  phosphorylated  in  fetal  progenitors  without  concurrent  MAPK  pathway 

262  activation.  This  raises  the  question  of  whether  each  pathway  has  unique 

263  functions  downstream  of  FLT3itd. 

264 

265  MAPK  pathway  inhibition  has  little  to  no  effect  on  Flt3ITD -driven  HSC  depletion, 

266  FI  PC  expansion  and  GMP  expansion. 

267  We  used  the  MEK  inhibitor  PD0325901  to  test  whether  MAPK  pathway 

268  inhibition  could  prevent  HSC  depletion  and  HPC/GMP  expansion  in  Flt3ITD  mice. 

269  We  administered  vehicle  or  PD0325901  to  6-week-old  wild  type  and  Flt3ITD/+ 

270  mice  (5  mg/kg  per  day  for  10  days).  This  regimen  effectively  inhibited  ERK1/2 

271  phosphorylation  in  HPCs  without  affecting  STAT5  phosphorylation  (Fig.  5  -  figure 

272  supplement  1).  PD0325901 -treated  wild  type  mice  had  significantly  more 

273  phenotypic  HSCs  and  HPCs  than  vehicle  treated  controls  (Fig.  5A,  B).  However, 

274  PD0325901  had  no  effect  on  HSC  numbers,  HPC  numbers  or  GMP  frequencies 
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275  in  Flt3ITD/+  mice  (Fig.  5A-C).  This  suggests  that  sustained  MAPK  pathway 

276  signaling  is  not  required  for  HSC  depletion,  HPC  expansion  and  GMP  expansion 

277  in  Flt3ITD/+  adult  mice. 

278  We  next  tested  whether  PD0325901  could  prevent  the  onset  of  the  HSC 

279  depletion,  HPC  expansion  and  GMP  expansion  phenotypes  if  it  was  given  shortly 

280  after  birth.  We  treated  nursing  mothers  of  wild  type,  Flt3ITD/+  and  Flt3ITD/ITD 

281  neonates  with  PD0325901  (5  mg/kg  per  day)  beginning  at  PI .  While  this  regimen 

282  has  previously  been  shown  to  rescue  MAPK  pathway-dependent  developmental 

283  abnormalities  in  Nfl  mutant  neonates  (Wang  et  al.,  2012),  it  did  not  prevent  HSC 

284  depletion  or  HPC  expansion  in  Flt3ITD  mutant  neonates  (Fig.  5D,  E),  and  it  only 

285  partially  rescued  GMP  expansion  (Fig.  5F).  Altogether,  the  data  suggest  that  the 

286  MAPK  pathway  has  only  a  minor  role,  if  any,  in  causing  these  phenotypes. 

287  Temporal  changes  in  MAPK  pathway  regulation  are  unlikely  to  account  for  the 

288  different  effects  of  FLT3itd  on  fetal  and  adult  progenitors. 

289 

290  STAT5  inactivation  exacerbates  FISC  depletion,  FiPC  expansion,  GMP 

291  expansion  and  MPN. 

292  STAT5  has  been  implicated  as  a  key  downstream  effector  of  FLT3itd  in 

293  many  different  systems,  and  it  is  hyper-phosphorylated  in  Flt3ITD/+  HSCs  and 

294  HPCs  during  fetal,  neonatal  and  adult  stages  of  development  (Fig.  4E,  F).  This 

295  raised  the  question  of  whether  genetic  inactivation  of  Stat5a  and  Stat5b  -  with  a 

296  conditional  Stat5a/b  allele  (Wang  et  al.,  2009)  -  could  prevent  HSC  depletion, 

297  HPC  expansion,  GMP  expansion  and  MPN  in  Flt3ITD/+  mice.  To  answer  this 
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298  question,  we  evaluated  HSCs,  HPCs,  GMPs  and  spleen  weights  in  1)  Flt3+/+ ; 

299  Stat5a/bf/+  or  Stat5a/bf/f ;  Cre  (control),  2)  Flt3+/+;  Stat5a/bf/+;  Mxl-Cre  ( Stat5 A/+), 

300  3)  Flt3+/+;  Stat5a/bf/f;  Mxl-Cre  ( Stat5 A/A),  4)  Flt3ITD/+ ;  Stat5a/bf/+  or  Stat5a/b1/f\  Cre 

301  ( Flt3ITD/+ ),  5)  Flt3ITD/+;  Stat5a/bf/f;  Mxl-Cre  ( Flt3ITD ;  Stat5 A/+),  and  6)  Flt3ITD/+; 

302  Stat5a/bf/f;  Mxl-Cre  ( Flt3ITD/+ ;  Stat5A/A)  mice.  The  mice  were  treated  with  poly- 

303  inosine:poly-cytosine  (plpC)  beginning  at  6  weeks  after  birth  to  delete  Stat5a/b, 

304  and  they  were  analyzed  4  weeks  later.  Western  blotting  confirmed  complete  loss 

305  of  STAT5  protein,  and  MAPK  pathway  activation  was  unaffected  by  Stat5a/b 

306  deletion  (Fig.  6  -  figure  supplement  1 ).  Surprisingly,  Stat5a/b  deletion 

307  exacerbated  the  HSC  depletion,  HPC  expansion  and  GMP  expansion 

308  phenotypes  of  Flt3ITD/+  mice  rather  than  rescuing  them  (Fig.  6A-C).  Spleen 

309  weights  were  also  enlarged  in  Flt3ITD;  Stat5A/+  and  Flt3ITD/+;  Stat5A/A  mice  relative 

310  to  control,  Stat5 'A/+,  Stat5A/A  and  Flt3ITD/+  littermates  (Fig.  5D).  Similar  results  were 

311  observed  when  we  deleted  a  single  Stat5a/b  with  Vavl-cre.  Only  one  Stat5a/b 

312  allele  was  deleted  in  these  analyses  because  bi-allelic  deletion  impairs  fetal 

313  erythropoiesis  (Zhu  et  al.,  2008).  Nevertheless,  Flt3ITD/+;  Stat5a/bf/+;  Vavl-Cre 

314  mice  had  fewer  HSCs,  more  HPCs,  more  GMPs  and  larger  spleens  than  control 

315  or  Flt3ITD/+  littermates  at  8-10  weeks  after  birth  (Fig.  6E-H). 

316  The  data  reveal  an  unanticipated  function  for  STAT5  in  pre-leukemic, 

317  Flt3ITD- mutant  progenitors.  They  suggest  that  STAT5  helps  to  maintain  Flt3ITD- 

318  mutant  HSCs  in  an  uncommitted  state  and  that  it  antagonizes  Flt3ITD- driven 

319  expansion  of  more  committed  myeloid  progenitor  populations.  Thus,  FLT3itd  may 
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320  simultaneously  potentiate  self-renewal  and  myeloid  commitment  programs  via 

321  STAT5-dependent  and  STAT5-independent  pathways,  respectively  (Fig.  61). 

322 

323  FLT3itd  activates  self-renewal  programs  in  HPCs  via  STAT5,  and  it  activates 

324  commitment  programs  independently  of  STAT5  The  changes  in  HSC  and 

325  HPC  frequencies  in  Flt3ITD;  Stat5a/b  compound  mutant  mice  raise  the  question  of 

326  whether  FLT3itd  has  STAT5-dependent  and  STAT5-independent  effects  on  gene 

327  expression,  and  whether  transcriptional  changes  are  developmental  context- 

328  specific.  To  answer  these  questions  we  performed  two  independent  experiments 

329  to  characterize  global  changes  in  gene  expression  (Fig.  7A).  In  the  first 

330  experiment,  we  analyzed  gene  expression  in  wild  type  and  Flt3ITD/+  HSCs  and 

331  HPCs  at  E14.5,  P0,  P14  and  8-10  weeks  after  birth.  This  experiment  was  meant 

332  to  elucidate  changes  in  FLT3itd  target  genes  over  time.  In  the  second 

333  experiment,  we  analyzed  gene  expression  in  adult  HPCs  from  1)  wild  type,  2) 

334  Flt3ITD/+,  3)  Flt3ITD/ITD,  4)  Flt3ITD;  Stat5A/+  and  5)  Flt3ITD/+;  Stat5A/A  mice.  This 

335  experiment  was  meant  to  delineate  which  FLT3itd  targets  are  STAT5-dependent 

336  and  which  are  STAT5-independent. 

337  We  analyzed  the  data  from  each  experiment  independently,  and  we 

338  merged  the  data  to  identify  a  list  of  genes  that  were  significantly,  differentially 

339  expressed  in  FLT3itd  progenitors  in  both  experiments.  In  experiment  1,  we 

340  identified  254  annotated  coding  genes  that  were  differentially  expressed  between 

341  wild  type  and  Flt3ITD/+  HPCs  at  one  or  more  time  points  (Fig.  7  -  source  data 

342  table  1 ;  adjusted  p<0.05;  fold  change  >2).  We  did  not  identify  any  genes  that  met 
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343  these  stringent  filtering  criteria  in  HSCs,  though  statistically  significant  changes  in 

344  gene  expression  were  observed  when  specific  target  genes  (from  the  HPC  list) 

345  were  individually  interrogated  (p<0.05,  Fig.  7  -  source  data  table  1 ).  The 

346  differences  between  HSCs  and  HPCs  may  simply  reflect  differences  in  Flt3 

347  expression  (Fig.  1A),  though  it  is  also  possible  that  HSCs  with  the  strongest 

348  transcriptional  responses  to  FLT3itd  were  not  captured  in  our  microarray  assays 

349  because  they  differentiated.  Of  the  254  genes  that  were  differentially  expressed 

350  in  experiment  1 ,  58  unique  genes  were  also  differentially  expressed  between  wild 

351  type  and  Flt3ITD/+  HPCs  in  experiment  2  (Fig.  7B  and  Fig.  7  -  figure  supplement 

352  1).  Thirty-three  genes  were  expressed  at  higher  levels  in  Flt3ITD  HPCs  relative  to 

353  wild  type  HPCs,  and  25  genes  were  expressed  at  lower  levels  (Fig.  7B).  Of 

354  these,  35  normalized  when  Stat5a/b  was  deleted,  but  23  did  not.  FLT3itd 

355  therefore  has  both  STAT5-dependent  and  STAT5-independent  effects  on  gene 

356  expression,  and  these  effects  are  more  pronounced  in  HPCs  as  compared  to 

357  HSCs. 

358  We  tested  whether  FLT3itd  activates  self-renewal-  and  commitment- 

359  related  transcriptional  programs  via  STAT5-dependent  and  STAT5-independent 

360  mechanisms,  respectively,  as  predicted  by  our  phenotypic  assays  (Fig.  6).  We 

361  generated  self-renewal-related  and  commitment-related  gene  sets  by  comparing 

362  wild  type  HSCs  and  HPCs  using  the  data  collected  in  experiment  1  (Figure  7  - 

363  source  data  table  2).  We  then  used  Gene  Set  Enrichment  Analysis  (GSEA)  to 

364  compare  wild  type,  Flt3ITD/+  and  Flt3ITD/+;  Stat5A/A  HPCs  using  data  collected  in 

365  experiment  2  (Subramanian  et  al.,  2005).  Self-renewal-related  genes  were 
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366  enriched  in  Flt3ITD/+  HPCs,  and  commitment-related  genes  were  enriched  in  wild 

367  type  HPCs  (Fig.  7C).  This  suggests  that  the  FLT3itd  protein  can  activate  self- 

368  renewal  mechanisms  in  otherwise  non-self-renewing  HPCs.  Remarkably,  these 

369  effects  were  strongly  reversed  when  Stat5a/b  was  deleted  (Fig.  7C).  A  separately 

370  curated  self-renewal  gene  set  from  Ivanova  et  al.  was  similarly  enriched  in  wild 

371  type  and  Flt3ITD/+  HPCs  as  compared  to  Flt3ITD/+;  Stat5 HPCs  (Fig.  7D) 

372  (Ivanova  et  al.,  2002).  These  findings  are  consistent  with  a  model  in  which 

373  FLT3itd  signals  via  STAT5  to  ectopically  activate  self-renewal  programs  in  HPCs, 

374  and  it  simultaneously  promotes  myeloid  commitment  via  STAT5-independent 

375  mechanisms. 

376  To  better  understand  the  STAT5-independent  mechanisms  that  promote 

377  myeloid  commitment,  we  performed  GSEA  on  Flt3ITD/+  and  Flt3ITD/+;  Stat5A/A 

378  HPCs  with  curated  gene  sets  in  the  MSigDB  database  (Subramanian  et  al., 

379  2005).  The  most  significantly  enriched  gene  sets  in  Flt3ITD/+;  Stat5A/A  HPCs  were 

380  generally  associated  with  increased  inflammatory  cytokine  signaling  (Fig.  7E). 

381  This  finding  is  intriguing  because  several  prior  studies  have  linked  inflammatory 

382  cytokine  signaling  to  loss  of  adult  HSC  self-renewal  capacity  and  myeloid 

383  differentiation  (Baldridge  et  al.,  2010;  Essers  et  al.,  2009;  Pietras  et  al.,  2016).  Of 

384  note,  we  did  not  observe  changes  in  STAT 1 ,  STAT3  or  AKT  phosphorylation  in 

385  Flt3ITD/+;  Stat5A/A  HPCs  (Fig.  6  -  figure  supplement  1  and  data  not  shown). 

386  Additional  studies  are  still  needed  to  identify  the  signal  transduction  molecules 

387  that  mediate  FLT3ITD-driven  myeloid  commitment,  and  to  test  whether  changes  in 
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388  cytokine-related  gene  expression  are  a  cause  or  a  consequence  of  differentiation 

389  in  Flt3ITD/+;  Stat5 HPCs. 

390 

391  FLT3ITD-mediated  changes  in  gene  expression  correlate  temporally  with  the 

392  normal  transition  from  fetal  to  adult  transcriptional  programs 

393  FLT3itd  has  the  capacity  to  activate  functionally  relevant  signal 

394  transduction  pathways,  such  as  STAT5,  during  both  pre-  and  post-natal  stages  of 

395  development  (Fig.  4),  yet  HSC  and  HPC  phenotypes  were  only  observed  after 

396  birth.  This  raises  the  question  of  whether  pre-  and  post-natal  progenitors  have 

397  distinct  transcriptional  responses  to  FLT3itd.  We  evaluated  FLT3itd  target  gene 

398  expression  in  HPCs  at  E14.5,  P0,  P14  and  adulthood  (Fig.  8A).  We  found  that 

399  differences  between  wild  type  and  Flt3ITD/+  HPCs  were  more  evident  at  PI 4  and 

400  adulthood  than  at  E14.5  or  P0  (Fig.  8A  and  Fig.  8  -  figure  supplement  1 ).  This 

401  was  true  for  both  STAT5-dependent  targets,  e.g.  Socs2,  and  STAT5-independent 

402  targets,  e.g.  Ctsg  (Fig.  8B).  Thus,  fetal,  neonatal  and  adult  hematopoietic 

403  progenitors  have  distinct  transcriptional  responses  to  FLT3itd  signaling. 

404  To  better  understand  when  HSCs  and  HPCs  transition  from  fetal  to  adult 

405  transcriptional  programs,  we  analyzed  gene  expression  in  wild  type  cells  from 

406  experiment  1 .  We  identified  2627  unique  genes  (from  3005  different  probes)  that 

407  exhibited  significant  changes  in  gene  expression  in  HSCs  between  E14.5  and 

408  adulthood.  Of  the  228  most  differentially  expressed  transcripts  (from  the  top  250 

409  probes),  all  followed  a  consistent  trend  toward  increasing  (109  genes)  or 

410  decreasing  (119  genes)  expression  with  increasing  age,  and  most  were 
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411  differentially  expressed  in  both  HSCs  and  HPCs  (Fig.  8C,  Fig.  8  -  source  data 

412  table  1 ).  Among  these  genes  were  several  that  encode  transcription  factors  and 

413  RNA  binding  proteins  that  are  known  to  regulate  HSC  self-renewal,  including 

414  Lin28b,  Esrl,  Hmga2  and  Egrl  (Fig.  8C,  D).  Principal  component  analysis  and 

415  Euclidean  distance  measurements  showed  that  PI 4  HSCs  more  closely 

416  resembled  adult  HSCs  than  fetal  HSCs,  and  PO  HSCs  more  closely  resembled 

417  fetal  HSCs  (Fig.  8E).  Similar  associations  were  observed  for  HPCs  (Fig.  8  - 

418  figure  supplement  2).  The  data  show  that  HSCs  and  HPCs  begin  transitioning 

419  from  fetal  to  adult  transcriptional  programs  by  PI 4,  even  before  they  achieve 

420  quiescence  (Bowie  et  al.,  2006).  Furthermore,  the  data  suggest  that  HSCs  and 

421  HPCs  become  competent  to  express  (or  repress)  FLT3itd  target  genes  as  they 

422  transition  from  fetal  to  adult  transcriptional  programs. 

423 

424  Flt3ITD  and  Runxl  heterozygous  mutations  collaboratively  induce  changes  in 

425  gene  expression  in  a  developmental  context-dependent  manner 

426  Flt3ITD  and  Tet2  loss  of  function  mutations  have  recently  been  shown  to 

427  cooperatively  induce  changes  in  gene  expression  and  DNA  methylation  in  adult 

428  HPCs  that  are  not  observed  with  either  mutation  alone  (Shih  et  al.,  2015).  We 

429  tested  whether  Flt3ITD  and  Runxl  mutations  have  similar  cooperative  effects  on 

430  transcription  and  whether  the  effects  are  age-specific.  We  evaluated  gene 

431  expression  in  1 )  wild  type,  2)  Runxl N\  3)  Flt3ITD/+  and  4)  Flt3ITD/+;  Runx1N+  HPCs 

432  at  P0  and  PI 4.  At  PI 4,  we  identified  191  genes  that  were  significantly 

433  differentially  expressed  between  wild  type  and  Flt3ITD/+;  Runxl HPCs  (adj. 
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434  p<0.05;  fold  change  >3).  At  PO  only  8  genes  met  these  criteria,  7  of  which 

435  overlapped  with  the  PI 4  gene  list  (Fig.  9A  and  Fig.  9  -  source  data  table  1 ). 

436  GSEA  showed  significant  overlap  between  genes  that  were  differentially 

437  expressed  in  Flt3ITD/+;  Runx1A/+  HPCs  and  those  that  Shih  et  al.  found  to  be 

438  differentially  expressed  in  Flt3ITD;  Tet 2A/A  HPCs  (Shih  et  al.,  2015)  (Fig.  9B).  This 

439  suggests  that  FLT3itd  can  cooperate  with  diverse  mutations  to  induce  a 

440  conserved  set  of  target  genes. 

441  We  used  hierarchical  clustering  to  better  visualize  differences  in  gene 

442  expression  at  PO  and  PI 4.  This  approach  did  show  some  differences  between 

443  wild  type  and  Flt3ITD/+;  Runx1A/+  HPCs  at  PO  (Fig.  9C),  but  much  greater 

444  differences  were  observed  at  PI 4  for  most  target  genes  (Fig.  9C,  D).  Of  the  25 

445  most  differentially  expressed  genes  in  P14  Flt3ITD/+;  Runx1A/+  HPCs,  only  three  - 

446  Nov,  Bhlhe40  and  Mboat2  -  were  induced  equally  at  both  PO  and  PI 4  (Figs.  9D, 

447  Fig.  9  -  source  data  table  1).  The  remaining  genes  showed  only  a  partial  change 

448  in  expression  at  PO  (e.g.  Socs2 )  or  no  change  in  expression  (e.g.  Adgrel,  Dusp6, 

449  Gem,  Postn)  (Fig.  9D).  GSEA  also  showed  differences  in  the  transcriptional 

450  programs  of  PO  and  P14  Flt3ITD/+;  Runx1A/+  HPCs.  Several  Gene  Ontology  and 

451  Oncogenic  Signatures  gene  sets  were  enriched  in  P14  Flt3ITD/+;  Runx1A/+  HPCs 

452  relative  to  wild  type  controls,  the  most  significant  of  which  included  genes  that 

453  are  negatively  regulated  by  mTOR  (Majumder  et  al.,  2004),  genes  that  inhibit 

454  apoptosis  (Gene  Ontology)  and  c-myc  targets  (Bild  et  al.,  2006)  (Fig.  9E).  These 

455  gene  sets  were  not  significantly  enriched  in  P0  HSCs  with  the  exception  of  the 

456  anti-apoptotic  gene  set,  which  was  paradoxically  enriched  in  wild  type  HPCs 


29 


457  relative  to  Flt3ITD/+;  Runx1A/+  (Fig.  9E).  Altogether,  the  data  show  that  cooperating 

458  Flt3ITD  and  Runxl  mutations  -  and  likely  other  cooperating  mutations  -  have 

459  developmental  context-specific  effects  on  gene  expression. 

460 

461  Discussion 

462  Our  data  offer  a  potential  explanation  for  why  FLT3itd  mutations  are  rare 

463  in  infants  and  young  children  with  AML.  The  mutant  FLT3itd  protein  hyper- 

464  activates  STAT5  in  hematopoietic  progenitors  during  both  fetal  and  adult  stages 

465  of  development,  yet  its  effects  on  transcription  are  realized  selectively  in  adult 

466  progenitors.  Even  cooperative  interactions  between  Flt3ITD  and  Runxl 

467  heterozygous  mutations  are  blunted  during  fetal  and  perinatal  stages.  Thus, 

468  FLT3itd  mutations  may  occur  disproportionately  in  older  children  and  adults  with 

469  AML  because  they  are  less  able  to  activate  key  effectors  of  leukemogenesis 

470  during  earlier  stages  of  life. 

471 

472  Flt3ITD  has  distinct  STA  T5-dependent  and  STA  T5-independent  effects  on  self- 

473  renewal  and  myeloid  commitment,  respectively 

474  In  the  course  of  these  studies,  we  discovered  that  Stat5a/b  deletion 

475  exacerbates,  rather  than  rescues,  the  HSC  depletion,  HPC/GMP  expansion  and 

476  MPN  phenotypes  of  Flt3ITD/+  mice.  Furthermore,  FLT3itd  ectopically  activated 

477  STAT5-dependent  self-renewal  programs  in  HPCs  (Fig.  7).  These  unanticipated 

478  findings  suggest  that  FLT3itd  activates  both  STAT5-dependent  and  STAT5- 

479  independent  signal  transduction  pathways  and  that  these  pathways  have 
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480  opposing  effects  on  self-renewal  and  myeloid  commitment.  In  this  model,  STAT5- 

481  independent  myeloid  commitment  programs  outweigh  STAT5-dependent  self- 

482  renewal  programs  in  Flt3ITD  mutant  bone  marrow  so  the  HSC  pool  becomes 

483  depleted  (Fig.  61  and  Fig.  7).  Stat5a/b  deletion  can  shift  the  balance  further  in 

484  favor  of  differentiation,  though  cooperating  mutations  may  ultimately  allow 

485  STAT5-dependent  self-renewal  programs  to  predominate  in  transformed  AML 

486  cells. 

487  The  STAT5-independent  pathways  that  antagonize  HSC  self-renewal  and 

488  promote  myeloid  progenitor  expansion  remain  unclear.  While  the  MAPK  pathway 

489  was  hyper-activated  in  postnatal  Flt3ITD/+  HSCs  and  HPCs,  MEK  inhibition  did  not 

490  prevent,  or  even  reduce,  HSC  depletion  or  HPC  expansion  in  these  mice  (Fig.  5). 

491  Other  candidate  pathways,  including  STAT3  and  PI3K,  were  not  activated  by 

492  FLT3itd  (Fig.  4A).  It  is  possible  that  low  levels  of  signal  transduction  via  these 

493  pathways  were  undetectable  by  Western  blot  but  nevertheless  functionally 

494  important.  It  is  also  possible  that  an  un-interrogated  pathway,  such  as  NF-kB  or 

495  CDK1  (Gerloff  et  al.,  2015;  Radomska  et  al.,  2012),  could  promote  myeloid 

496  commitment  and  antagonize  STAT5.  Our  GSEA  data  did  show  increased 

497  expression  of  inflammatory  cytokine  receptors  in  Stat5a/b- deficient,  Flt3ITD  HPCs. 

498  This  raises  the  intriguing  possibility  that  inflammatory  cytokines  could  promote 

499  differentiation  of  Flt3ITD  mutant  progenitors,  and  perhaps  AML  cells.  Additional 

500  genetic  studies  are  needed  to  resolve  whether  these  transcriptional  changes  are 

501  a  cause  or  a  consequence  of  enhanced  lineage  commitment  in  Flt3ITD  mutant 

502  HPCs.. 
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503 

504  Developmental  programming,  re-programming  and  the  origins  of  pediatric  and 

505  adult  malignancies 

506  Heterochronic  genes  have  been  implicated  in  cancer  pathogenesis  (Shyh- 

507  Chang  and  Daley,  2013).  For  example,  hepatoblastomas,  Wilm’s  tumors  and 

508  most  neuroblastomas  present  early  in  life,  and  they  often  express  the  oncofetal 

509  proteins  LIN28  or  LIN28B  at  high  levels  to  help  maintain  the  primitive 

510  differentiation  states  of  their  respective  anlagen  (Diskin  et  al.,  2012;  Molenaar  et 

511  al.,  2012;  Nguyen  et  al.,  2014;  Urbach  et  al.,  2014).  Adult  hepatocellular 

512  carcinomas,  germ  cell  tumors  and  ovarian  carcinomas  often  ectopically  activate 

513  LIN28  or  LIN28B  to  restore  oncofetal  programs  (Viswanathan  et  al.,  2009).  MLL- 

514  rearranged  leukemias  have  similarly  been  shown  to  express  embryonic  stem 

515  cell-related  genes  (Somervaille  et  al.,  2009),  and  BRAFV600E  driven  melanoma 

516  was  recently  shown  to  arise  from  melanocytes  that  first  de-differentiate  into 

517  primitive  neural  crest  progenitors  (Kaufman  et  al.,  2016).  In  each  of  these  cases, 

518  it  is  easy  to  appreciate  why  maintaining  or  restoring  a  primitive  cell  identity  might 

519  accelerate  transformation  -  fetal  cells  can  proliferate  rapidly  without 

520  differentiating  or  senescing.  However,  it  is  then  curious  as  to  why  pediatric 

521  malignancies  are  relatively  uncommon.  Does  this  simply  reflect  greater  fidelity  of 

522  the  genome  at  early  ages,  or  are  other  factors  at  work? 

523  Our  data  suggest  that  the  transcriptional  regulatory  programs  of  fetal 

524  progenitors  may,  in  fact,  be  protective  against  some  mechanisms  of 

525  transformation.  Fetal  and  adult  progenitors  interpret  FLT3ITD-derived  signals 
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526  differently,  as  evidenced  by  their  distinct  transcriptional  responses  (Figs.  7  and 

527  8),  and  this  constrains  the  ability  of  FLT3itd  to  transform  fetal  progenitors.  Either 

528  they  lack  key  transcriptional  co-activators,  or  the  epigenetic  landscape  of  fetal 

529  progenitors  suppresses  FLT3itd  target  gene  activation.  Further  work  is  needed  to 

530  understand  the  cis-  and  trans-regulatory  elements  that  determine  when  and  how 

531  individual  mutations  are  competent  to  transform.  If  we  can  understand  how 

532  normal  developmental  programs  interact  with  genetic  mutations  to  cause 

533  malignancies,  it  may  be  possible  to  target  these  interactions  therapeutically. 

534 

535  Materials  and  methods 

536  Mouse  strains 

537  The  Flt3ITD  RRID:IMSR_JAX:01 1112  (Lee  et  al„  2007),  Runx1f 

538  RRID:IMSR_JAX:008772  (Taniuchi  et  al„  2002),  Stat5a/bf 

539  RRI D: M M RRC  032053- J  AX  (Cui  et  al„  2004),  Rictor1  RRID:IMSR_JAX:020649 

540  (Magee  et  al.,  2012),  Vavl-Cre  RRID:IMSR_JAX:008610  (Siegemund  et  al., 

541  2015)  and  Mxl-Cre  RRID:IMSR_JAX:003556  (Kuhn  et  al.,  1995)  mouse  strains 

542  have  all  been  previously  described  and  were  obtained  from  The  Jackson 

543  Laboratory.  These  lines  were  all  on  a  pure  C57BL/6  background.  Expression  of 

544  Mxl-Cre  was  induced  by  three  intraperitoneal  injections  of  pIpC  (GE  Life 

545  Sciences;  10  j^g/dose)  over  five  days  beginning  6  weeks  after  birth.  PD0325901 

546  (Cayman  Chemicals)  was  suspended  in  0.5%  hydroxypropylmethylcellulose 

547  vehicle  (Sigma)  and  administered  by  oral  gavage  as  described  in  the  text.  All 

548  mice  were  housed  in  the  Department  for  Comparative  Medicine  at  Washington 
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549  University.  All  animal  procedures  were  approved  by  the  Washington  University 

550  Committees  on  the  Use  and  Care  of  Animals. 

551 

552  Isolation  of  HSCs  and  flow  cytometry 

553  Bone  marrow  cells  were  obtained  by  flushing  the  long  bones  (tibias  and 

554  femurs)  or  by  crushing  long  bones,  pelvic  bones  and  vertebrae  with  a  mortar  and 

555  pestle  in  calcium  and  magnesium-free  Hank’s  buffered  salt  solution  (HBSS), 

556  supplemented  with  2%  heat  inactivated  bovine  serum  (Gibco).  Splenocytes  were 

557  obtained  by  macerating  spleens  with  frosted  slides.  Single  cell  suspensions  were 

558  filtered  through  a  40  um  cell  strainer  (Fisher).  The  cells  were  then  stained  for  20 

559  minutes  with  fluorescently  conjugated  antibodies,  washed  with  HBSS  +  2% 

560  bovine  serum  and  resuspended  for  analysis.  Cell  counts  were  measured  by 

561  hemocytometer.  The  following  antibodies  were  used  for  flow  cytometry,  all  were 

562  from  Biolegend  except  as  indicated:  CD150  (TC15-12F12.2),  CD48  (HM48-1), 

563  Seal  (D7),  c-Kit  (2B8),  Ter119  (Ter-119),  CD3  (17A2),  CDIIb  (Ml/70),  Gr-1 

564  (RB6-8C5),  B220  (RA3-6B2),  CD8a  (53-6.7),  CD34  (eBioscience,  RAM 34),  CD2 

565  (RM2-5),  CD45.1  (A20),  CD45.2  (104),  CD127  (A7R340),  CD16/32  (93)  and 

566  FLT3/CD135  (A2F10).  Lineage  stains  for  all  experiments  included  CD2,  CD3, 

567  CD8a,  Ter119,  B220  and  Grl.  Antibodies  to  CD4  and  CDIIb  were  omitted  from 

568  the  lineage  stains  because  they  are  expressed  on  fetal  HSCs  at  low  levels. 

569  Unless  otherwise  indicated,  we  used  the  following  surface  marker  phenotypes  to 

570  define  cell  populations:  HSCs  (CD150+,  CD48'Lineage',  Sca1+,  c-kit+),  HPCs 

571  (CD48'Lineage',  Sca1+,  c-kit+),  and  GMPs  (Lineage',  Seal',  CD127',  c-kit+, 
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572  CD34+,  CD16/32+).  Non-viable  cells  were  excluded  from  analyses  by  4’, 6- 

573  diamidino-2-phenylindone  (DAPI)  staining  (1  jig/ml).  When  HSCs  and  HPCs 

574  were  isolated  for  Western  blotting  or  RNA  collection,  c-kit+  cells  were  enriched 

575  prior  to  sorting  by  selection  with  paramagnetic  beads  (Miltenyi  Biotec).  Flow 

576  cytometry  was  performed  on  a  BD  FACSAria  Fusion  flow  cytometer  (BD 

577  Biosciences). 

578 

579  Limiting  dilution  long-term  repopulation  assays 

580  Eight  to  ten  week  old  C57BL/6Ka-Thy-1 .2  (CD45.1 )  recipient  mice  were 

581  given  two  doses  of  550  rad  delivered  at  least  3  hours  apart.  Donor  fetal  liver  or 

582  bone  marrow  cells  were  mixed  with  competitor  bone  marrow  cells  at  the  doses 

583  indicated  in  the  text  and  injected  via  the  retroorbital  sinus.  To  assess  donor 

584  chimerism,  peripheral  blood  was  obtained  from  the  submandibular  veins  of 

585  recipient  mice  at  the  indicated  times  after  transplantation.  Blood  was  subjected  to 

586  ammonium-chloride  lysis  of  the  red  blood  cells  and  leukocytes  were  stained  with 

587  antibodies  to  CD45.2,  CD45.1 ,  B220,  CD3,  CD1 1  b  and  Gr-1  to  assess 

588  multilineage  engraftment.  Functional  HSC  frequencies  were  calculated  and 

589  compared  by  using  Extreme  Limiting  Dilution  Analysis  (Hu  and  Smyth,  2009).  For 

590  secondary  transplants,  mice  were  injected  with  3  million  cells  from  the  bone 

591  marrow  of  primary  recipient  mice. 

592 

593 

594  Quantitative  RT-PCR 

595  RNA  was  isolated  from  HSCs  with  RNAeasy  micro  plus  columns  (Qiagen) 

596  and  converted  to  cDNA  with  Superscript  III  reverse  transcriptase  (Lifetech). 
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597  Quantitative  RT-PCR  assays  were  performed  with  Taqman  Gene  Expression 

598  Assays  specific  to  mouse  Flt3  and  p-actin  (Lifetech).  Analysis  was  performed 

599  with  a  Mx3005P  qPCR  system  (Agilent).  Samples  were  normalized  based  on  |3- 

600  actin  expression. 

601 

602  Western  blots 

603  Twenty-five  thousand  HSC/MPPs,  HPCs  or  GMPs  were  double  sorted  into 

604  Trichloracetic  acid  (TCA),  and  the  volume  was  adjusted  to  a  final  concentration  of 

605  10%  TCA.  Extracts  were  incubated  for  15  minutes  on  ice  and  centrifuged  at 

606  16,100xg  at  4°C  for  10  minutes.  Precipitates  were  washed  in  acetone  twice  and 

607  dried.  The  pellets  were  solubilized  in  9M  urea,  2%  Triton  X-100,  1%  DTT.  LDS 

608  loading  buffer  (Lifetech)  was  added  and  the  pellet  was  heated  at  70°C  for  10 

609  minutes.  Samples  were  separated  on  Bis-Tris  polyacrylamide  gels  (Lifetech)  and 

610  transferred  to  PVDF  membrane  (Lifetech).  All  antibodies  were  from  Cell 

611  Signaling  Technologies  except  as  indicated:  P-STAT5  (4322),  Total  STAT5 

612  (9363),  P-STAT3  (9145),  Total  STAT3  (9139),  P-ERK1/2  (4370),  Total  ERK1/2 

613  (4696),  P-AKT  Ser473  (4060),  P-AKT  T308  (13038),  Total  AKT  (4691 ),  cx- 

614  TUBULIN  (3873),  p-ACTIN  (Santa  Cruz  Bioscience,  clone  AC-17),  HRP-anti- 

615  Rabbit  IgG  (7074)  and  HRP-anti-mouse  IgG  (7076).  Blots  were  developed  with 

616  the  SuperSignal  West  Femto  or  Pico  chemiluminescence  kits  (Thermo  Scientific). 

617  Blots  were  stripped  (1%  SDS,  25  mM  glycine  pH  2)  prior  to  re-probing. 

618 

619  Cytospins 
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620  Bone  marrow  cells  were  isolated  and  spun  onto  glass  slides  using  a 

621  Shandon  Cytospin  3.  The  slides  were  stained  using  Protocol  Hema  3  Wright- 

622  Giemsa  stain  (Fisher  Scientific).  All  slides  were  reviewed  by  a  pediatric 

623  hematologist  (JAM  or  ASC). 

624 

625  Statistical  analysis 

626  In  all  cases,  multiple  independent  experiments  were  performed  on  at  least 

627  three  different  days  to  verify  that  the  data  are  reproducible.  Grouped  data  reflect 

628  biological  replicates  (i.e.  independent  mice)  and  are  represented  by  mean  +/- 

629  standard  deviation.  Statistical  comparisons  between  groups  were  made  with  the 

630  two-tailed  Student’s  t-test  except  as  noted  in  the  figure  legends.  When  multiple 

631  genotypes  were  compared,  statistical  significance  was  determined  by  performing 

632  a  one-way  ANOVA  followed  by  a  Holm-Sidak  post-hoc  test  to  correct  for  multiple 

633  comparisons.  For  transplantation  experiments,  the  percentages  of  mice  with 

634  multilineage  reconstitution  were  compared  with  the  Fisher’s  exact  test.  All 

635  comparisons  were  performed  with  GraphPad  Prism  6,  RRID:SCR_002798. 

636 

637  Gene  expression  analysis 

638  Ten  thousand  HSCs  or  HPCs  were  double  sorted  directly  into  RLT  plus 

639  RNA  lysis  buffer  (Qiagen)  and  RNA  was  isolated  with  RNAeasy  micro  plus 

640  columns  (Qiagen).  Transcripts  were  amplified  with  the  WTA2  kit  (Sigma)  with  the 

641  Kreatech  ULS  RNA  labeling  kit  (Kreatech  Diagnostics).  Labeled  cDNA  was 

642  hybridized  to  Agilent  Mouse  8x60K  microarrays  and  analyzed  with  an  Agilent  C- 
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643  class  scanner.  Signal  data  were  assembled  and  processed  in  Partek 

644  RRID:SCR_01 1860,  and  samples  were  compared  by  Linear  Models  for 

645  Microarrays,  RRID:SCR_010943  (Ritchie  et  al.,  2015;  Smyth,  2004).  Adjusted  p- 

646  values  were  calculated  by  the  Benjamini  and  Hochberg  false  discovery  rate 

647  (Benjamini  and  Hochberg,  1995).  Z-scores  were  calculated  as  previously 

648  described  (Cheadle  et  al.,  2003).  Hierarchical  cluster  analysis  was  performed 

649  with  Cluster  3.0  and  visualized  with  Java  TreeView;  RRI D : SC R_0 13505  and 

650  RRID:SCR_013503  (Eisen  et  al.,  1998).  Principal  component  analyses  and 

651  Euclidean  distance  comparisons  (by  permutation  testing)  were  performed  with 

652  the  R  software  environment.  Microarray  data  sets  have  been  deposited  into 

653  Gene  Expression  Omnibus  (GSE81 153).  GSEA  was  performed  using  gene  sets 

654  that  were  generated  as  cited  in  the  text,  or  with  gene  sets  curated  in  the  MSigDB 

655  databases;  RRID:SCR_003199  (Subramanian  et  al.,  2005). 

656 
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916  Figure  Legends 

917 

918  Figure  1.  Flt3ITD  causes  HSC  depletion  in  adult  but  not  fetal  mice. 

919  (A)  Flt3  transcript  expression  in  fetal  and  adult  HSCs  and  HPCs  relative  to  fetal 

920  HSCs;  n=4-9.  (B)  FLT3  expression  in  fetal  and  adult  HSC/HPCs  (Lineage' 

921  Sca1+c-kit+)  and  unfractionated  fetal  liver  or  bone  marrow  cells,  as  determined  by 

922  flow  cytometry  (N=3).  (C)  HSC  numbers  in  two  tibias  and  femurs  from  adult  wild 

923  type  and  Flt3ITD  mice;  n=12-16.  (D)  HSC  numbers  in  fetal  livers  from  El 4.5  wild 

924  type  and  Flt3ITD  mice;  n=9-20.  (E)  Spleen  HSC  frequency  in  adult  wild  type  and 

925  Flt3ITD/ITD  mice;  n=4-5.  (F,  G)  Limiting  dilution  analyses  using  adult  bone  marrow 

926  (F)  or  E14.5  fetal  liver  cells  (G);  n=9-10  recipients  per  cell  dose.  Wild  type  and 

927  Flt3ITD/ITD  HSC  frequencies  were  calculated  by  extreme  limiting  dilution  analysis. 

928  (H,  I)  Frequencies  of  donor  (CD45.2)  and  competitor  (CD45.1 )  HSCs  (H)  and 

929  donor  bone  marrow  cells  (I)  in  primary  recipients  of  100,000  fetal  liver  cells;  n=15 

930  per  genotype.  (J)  Frequencies  of  CD45.2+  peripheral  blood  cells  in  secondary 

931  recipients  of  donor  cells  that  originated  from  wild  type  or  Flt3ITD/ITD fetal  livers; 

932  n=12-14.  (K)  Percentage  of  secondary  recipient  mice  with  multilineage  donor 

933  reconstitution.  In  all  panels,  error  bars  indicate  standard  deviations  and  n  reflects 

934  biological  replicates.  *  p<0.05,  ***  p<0.001  by  two-tailed  Student’s  t-test.  ## 

935  p<0.01  by  Fisher  exact  probability  test. 

936 

937  Figure  2.  Flt3ITD  causes  HSC  depletion,  HPC  expansion  and  GMP  expansion 

938  at,  or  shortly  after,  birth. 
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939  (A)  Absolute  HSC  numbers  in  fetal  or  PO  livers  for  the  indicated  genotypes.  (B) 

940  Absolute  HSC  numbers  in  P14  and  adult  bone  marrow  (2  hind  limbs)  or  P14 

941  spleen.  (C,  D)  Fetal  and  adult  HPC  numbers  (2  hind  limbs).  (E)  GMP  frequencies 

942  in  fetal  liver  or  adult  bone  marrow.  (F)  Liver  or  spleen  weights.  In  all  panels,  error 

943  bars  indicate  standard  deviations;  n=6-20  biological  replicates  for  each  age  and 

944  genotype.  *  p<0.05;  **  p<0.01;  ***  p<0.001  by  two-tailed  Student’s  t-test  relative 

945  to  the  wild  type  control  at  the  same  time  point. 

946 

947  Figure  3.  Flt3ITD  and  Runxl  mutations  cooperate  to  deplete  HSCs  and 

948  expand  committed  progenitor  populations  after  birth. 

949  (A-D)  HSC  frequencies  in  E14.5  fetal  liver,  PO  liver,  P14  bone  marrow  and  P21 

950  bone  marrow  for  the  indicated  genotypes.  (E-H)  HPC  frequencies  in  E14.5  fetal 

951  liver,  PO  liver,  P14  bone  marrow  and  P21  bone  marrow  for  the  indicated 

952  genotypes.  (I-L)  GMP  frequencies  in  E14.5  fetal  liver,  PO  liver,  P14  bone  marrow 

953  and  P21  bone  marrow  for  the  indicated  genotypes.  (M,  N)  Percentages  of 

954  CD45.2+  donor  leukocytes  (M)  or  CD1 1  b+Gr1 +  myeloid  cells  (N)  in  the  peripheral 

955  blood  of  recipients  of  PO  liver  or  P21  bone  marrow  cells  from  control  or  Flt3ITD/+: 

956  Runx1N+  mice.  Measurements  are  shown  at  2  and  4  weeks  after  transplantation. 

957  (O,  P)  Percentage  of  recipients  with  multi-lineage  (O)  or  myeloid  (P)  donor 

958  reconstitution  at  4  weeks  after  transplantation.  In  all  panels,  error  bars  indicate 

959  standard  deviations.  For  A-L,  n=8-18  biological  replicates  per  genotype  and  age. 

960  For  M-P,  n=14-15  recipients  from  3  independent  donors.  Statistical  significance 

961  was  determined  with  a  one-way  ANOVA  followed  by  Holm-Sidak’s  post-hoc  test 
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962  for  multiple  comparisons  (*  p<0.05;  **  p<0.01 ;  ***  p<0.001 ),  or  #  p<0.0001  by  the 

963  Fisher  exact  probability  test. 

964 

965  Figure  4.  FLT3itd  activates  STAT5  in  both  fetal  and  adult  progenitors,  but  it 

966  activates  the  MAPK  pathway  after  birth. 

967  (A)  Western  blot  showing  phosphorylation  of  STAT5,  ERK1/2,  STAT3  and  AKT  in 

968  adult  wild  type  and  Flt3ITD/+  HSC/MPPs,  HPCs  and  GMPs.  (B)  Flt3ITD/+;  Runx1A/A 

969  progenitors  give  rise  to  AML  in  adult  mice  (right  panel)  that  is  not  observed  in 

970  wild  type  or  Flt3ITD/+  bone  marrow.  Scale  bars  indicate  100  microns.  (C)  STAT5 

971  and  MAPK  are  hyper-phosphorylated  in  Flt3ITD/+;  Runx1NA  AML  that  develops  in 

972  adult  mice.  (D)  STAT5  and  ERK1/2  phosphorylation  in  wild  type  and  Flt3ITD/+ 

973  HSC/MPPs  at  E14.5,  P0,  P14  and  adulthood.  (E)  STAT5  phosphorylation  in  wild 

974  type  and  Flt3ITD/+  HPCs  at  E14.5,  P0,  P14  and  adulthood.  (F)  ERK1/2 

975  phosphorylation  in  wild  type  and  Flt3ITD/+  HPCs  at  E14.5,  P0,  P14  and  adulthood. 

976  Each  blot  is  representative  of  2  (panels  A  and  C)  or  at  least  3  (panels  D-F) 

977  independent  experiments. 

978 

979  Figure  5.  MAPK  pathway  inhibition  does  not  prevent  HSC  depletion  or 

980  committed  progenitor  expansion  in  Flt3ITD/+  mice. 

981  (A-D)  HSC  numbers  (A),  HPC  numbers  (B)  and  GMP  frequencies  (C)  in  wild  type 

982  and  Flt3ITD/+  mice  that  were  treated  with  vehicle  or  PD0325901  for  10  days 

983  beginning  at  6  weeks  after  birth;  n=4-5  biological  replicates  per  genotype  and 

984  treatment.  (D-F)  HSC  numbers  (D),  HPC  numbers  (E)  and  GMP  frequencies  (F) 
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985  in  PI 9  wild  type,  Flt3ITD/+  and  Flt3ITD/ITD  mice  whose  mothers  were  given 

986  PD0325901  beginning  at  PI ;  n=4-1 5  biological  replicates  for  each  genotype  and 

987  treatment.  In  all  panels,  error  bars  indicate  standard  deviation.  Statistical 

988  comparisons  were  made  with  a  two-tailed  Student’s  t-test.  *  p<0.05  relative  to 

989  vehicle  treated  cells  with  equivalent  genotypes;  #  p<0.05  relative  to  similarly 

990  treated  wild  type  controls;  A  p<0.05  relative  to  similarly  treated  wild  type  and 

991  Flt3ITD/+  groups. 

992 

993  Figure  6.  Stat5a/b  deletion  exacerbates  rather  than  rescues  HSC  depletion, 

994  HPC  expansion,  GMP  expansion  and  MPN  in  F/f3,ro/+mice. 

995  (A-D)  HSC  numbers  (A),  HPC  numbers  (B),  GMP  frequencies  (C)  and  spleen 

996  weights  (D)  in  Flt3ITD/+;  Stat5a/bf/f;  Mxl-Cre  compound  mutant  mice  and  littermate 

997  controls;  n=6-20  biological  replicates  per  genotype.  Stat5a/b  was  conditionally 

998  deleted  6  weeks  after  birth,  and  this  caused  a  complete  loss  of  protein 

999  expression  (figure  supplement  1).  (E-H)  HSC  numbers  (E),  HPC  numbers  (F), 

1000  GMP  frequencies  (G)  and  spleen  weights  (H)  in  Flt3ITD/+;  Stat5a/bf/+;  Vavl-Cre 

1001  compound  mutant  mice  and  littermate  controls;  n=8-20  biological  replicates  per 

1002  genotype.  (I)  The  data  suggest  that  STAT5-dependent  pathways  promote  HSC 

1003  self-renewal  downstream  of  FLT3itd,  but  these  effects  are  outweighed  by  STAT5- 

1004  independent  myeloid  commitment  pathways.  In  all  panels,  error  bars  indicate 

1005  standard  deviation.  Statistical  significance  was  determined  with  a  one-way 

1006  ANOVA  followed  by  Holm-Sidak’s  post-hoc  test  for  multiple  comparisons. 

1007  *p<0.05;  **  p<0.01 ;  ***  p<0.001 . 
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1008 


1009  Figure  7.  FLT3itd  activates  STAT5-dependent  self-renewal  programs  and 

1010  STAT5-independent  commitment  programs. 

1011  (A)  Overview  of  experimental  design.  (B)  Heatmap  representing  genes  that  were 

1012  differentially  expressed  in  Flt3ITD  mutant  HPCs  relative  to  wild  type  HPCs  in  both 

1013  experiments  1  and  2.  Each  column  represents  an  independent  sample.  The  gene 

1014  names  and  dendrogram  are  shown  in  figure  supplement  1  attached  to  this  figure. 

1015  (C)  Self-renewal  and  commitment-related  gene  sets  were  generated  by 

1016  identifying  genes  that  were  more  highly  expressed  (>5  fold,  adj.  p<0.05)  in  HSCs 

1017  relative  to  HPCs  (self-renewal),  or  HPCs  relative  to  HSCs  (commitment).  GSEA 

1018  plots  show  ectopic  activation  of  self-renewal-related  genes  in  HPCs  that  express 

1019  FLT3itd,  but  these  effects  are  reversed  in  Stat5a/b- deficient  HPCs.  (D)  An 

1020  independently  curated  self-renewal  gene  set  (Ivanova  et  al.,  2002)  was  similarly 

1021  enriched  in  wild  type  and  Flt3ITD/+  HPCs  relative  to  Flt3ITD/+;  Stat5A/A  HPCs.  (E) 

1022  GSEA  revealed  enrichment  of  gene  sets  associated  with  increased  inflammatory 

1023  cytokine  signaling. 

1024 

1025  Figure  8.  F/f3/TD-mediated  changes  in  gene  expression  correlate  with  the 

1026  normal  transition  from  fetal  to  adult  transcriptional  programs. 

1027  (A)  Heatmap  showing  expression  of  FLT3itd  target  genes  at  E14.5,  P0,  P14  and 

1028  adult  stages.  Each  column  shows  average  fold  change  in  Flt3ITD/+  HPCs  relative 

1029  wild  type  HPCs  at  the  indicated  time  point;  n=3-4  independent  arrays  per 

1030  genotype.  The  gene  names  and  dendrogram  are  shown  in  figure  supplement  1 
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1031  attached  to  this  figure.  (B)  Representative  examples  of  expression  of  STAT5- 

1032  independent  ( Ctsg )  and  STAT5-dependent  ( Socs2 )  FLT3itd  targets.  Error  bars 

1033  reflect  standard  deviation.  ***adj.  p<0.05  relative  to  wild  type  at  the  same  time 

1034  point,  #  adj.  p<0.05  relative  to  Flt3ITD/+  at  the  same  time  point.  (C)  Heterochronic 

1035  genes  began  transitioning  from  fetal  to  adult  expression  patterns  between  P0  and 

1036  PI 4,  concordant  with  sensitivity  to  FLT3itd.  Genes  that  encode  transcription 

1037  factors  and  RNA  binding  proteins  are  noted  to  the  right  of  the  heatmap.  A 

1038  complete  gene  list  is  provided  in  source  data  table  1  attached  to  this  figure.  (D) 

1039  Representative  examples  of  heterochronic  genes  that  show  decreased  ( Lin28b ) 

1040  or  increased  (Esrl;  Estrogen  Receptor  a)  expression  in  adult  relative  to  fetal 

1041  HSCs  and  HPCs.  Error  bars  reflect  standard  deviation.  $$  adj.  p<0.05  relative  to 

1042  E14.5  for  both  HSCs  and  HPCs.  (E)  Principal  component  analysis  and  Euclidean 

1043  distance  measurements  show  that  gene  expression  in  P0  HSCs  more  closely 

1044  resembles  fetal  HSCs  than  adult  HSCs,  and  gene  expression  in  PI 4  HSCs  more 

1045  closely  resembles  that  of  adult  HSCs.  Similar  calculations  for  HPCs  are  shown  in 

1046  figure  supplement  2. 

1047 

1048  Figure  9.  Flt3ITD  and  Runxl  mutations  cooperatively  induce  changes  in 

1049  gene  expression  at  PI 4,  yet  they  have  a  much  smaller  effect  at  P0. 

1050  (A)  Venn  diagram  showing  overlap  between  genes  that  were  significantly 

1051  differentially  expressed  (adj.  p<0.05,  fold  change  >3)  in  Flt3ITD/+;  Runxl A/+  HPCs 

1052  relative  to  wild  type  at  PI 4  or  P0;  n=4-5  arrays  per  genotype  and  age.  (B)  GSEA 

1053  shows  that  differentially  expressed  genes  in  Flt3ITD/+;  Runx1A/+  HPCs  overlap 
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1054  significantly  with  genes  that  are  differentially  expressed  in  Flt3ITD/+;  Tet2A/A  HPCs 

1055  (Shih  et  al.,  201 5).  (C)  Heatmap  showing  expression  of  genes  that  were 

1056  differentially  expressed  in  Flt3ITD/+;  Runx1A/+  HPCs  relative  to  wild  type  HPCs. 

1057  Each  column  indicates  average  fold  change  relative  to  the  wild  type  samples 

1058  from  the  same  time  point.  The  gene  list  is  shown  in  source  data  table  1  attached 

1059  to  this  figure.  (D)  Representative  examples  of  genes  that  are  among  the  most 

1060  differentially  expressed  in  Flt3ITD/+;  Runx1A/+  HPCs  relative  to  wild  type  HPCs  at 

1061  P14.  Most  show  much  smaller  changes  in  expression  at  P0.  Error  bars  reflect 

1062  standard  deviations,  *  adj.  p<0.05.  (E)  GSEA  identified  several  gene  sets  that 

1063  were  enriched  in  Flt3ITD/+;  Runx1A/+  HPCs  relative  to  wild  type  HPCs  at  PI 4. 

1064  Three  of  the  most  significantly  enriched  gene  sets  are  shown  for  P14  and  P0. 

1065 

1066  Supplemental  Figure  Legends 

1067 

1068  Figure  4  -  figure  supplement  1.  The  PI3K/mTORC2  pathway  does  not 

1069  mediate  HSC  depletion  or  MPN  in  Flt3ITD  mutant  mice. 

1070  (A)  Western  blot  showing  STAT5  and  AKT  phosphorylation  in  wild  type  and 

1071  Flt3ITD/ITD  HPCs.  Pten- deficient  HPCs  were  included  as  a  positive  control  for  AKT 

1072  phosphorylation.  Rictor  deletion  prevented  AKT  hyper-phosphorylation  by 

1073  mTORC2.  (B,C)  HSC  numbers  in  two  tibias  and  femurs  (B)  and  spleen  weights 

1074  (C)  in  Flt3ITD/ITD;  Ricto //f;  Vavl-Cre  compound  mutant  mice  (and  the  indicated 

1075  controls)  after  conditional  Rictor  deletion;  n=3-12  per  genotype.  Error  bars 

1076  indicate  standard  deviations. 


59 


1077 


1078  Figure  5  -  figure  supplement  1 .  Inhibition  of  the  MAPK  pathway  fails  to 

1079  rescue  FLT3ITD-mediated  HSC  depletion  and  myeloid  progenitor  expansion, 

1080  but  Stat5a/b  deletion  enhances  these  phenotypes. 

1081  A  Western  blot  was  performed  with  25,000  HPCs  from  wild  type  and  Flt3ITD/+, 

1082  vehicle  and  PD0325901  treated  mice.  PD0325901  inhibited  ERK1/2 

1083  phosphorylation  without  affecting  STAT5  phosphorylation. 

1084 

1085  Figure  6  -  figure  supplement  1.  Stat5a/b  deletion  causes  a  complete  loss  of 

1086  phosphorylated  and  total  STAT5  protein. 

1087  A  Western  blot  was  performed  with  25,000  HPCs  from  wild  type,  Flt3ITD/+  and 

1088  Flt3ITD/+;  Stat5A/A  adult  mice  (representative  of  2  independent  experiments). 

1089  ERK1/2  phosphorylation  was  not  affected  by  STAT5  deletion.  STAT3 

1090  phosphorylation  was  not  affected,  and  STAT1  and  AKT  phosphorylation  were  not 

1091  detectable  (not  shown). 

1092 

1093  Figure  7  -  figure  supplement  1.  FLT3itd  induces  STAT5-dependent  and 

1094  STAT5-independent  changes  in  gene  expression. 

1095  An  expanded  version  of  Figure  7B  with  the  dendrogram  and  gene  names 

1096  attached  to  the  heatmap. 

1097 
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1098  Figure  8  -  figure  supplement  1.  FLT3ITD-mediated  changes  in  gene 

1099  expression  correlate  temporally  with  a  transition  from  fetal  to  adult 

1100  transcriptional  states. 

1101  An  expanded  version  of  Figure  8A  with  the  dendrogram  and  gene  names 

1102  attached  to  the  heatmap. 

1103 

1104  Figure  8  -  figure  supplement  2.  HPCs  express  heterochronic  genes  and 

1105  begin  to  transition  from  fetal  to  adult  transcriptional  programs  by  PI 4. 

1106  Principal  component  analysis  shows  temporal  changes  in  HPC  gene  expression 

1107  before  and  after  birth.  Euclidean  distance  measurements  and  permutation  testing 

1108  show  that  P0  HPCs  more  closely  resemble  fetal  HPCs  than  adult  HPCs.  P14 

1109  HPCs  more  closely  resemble  adult  HPCs. 

1110 
1111 

1112  Source  Data  Tables 

1113 

1114  Figure  7  -  source  data  table  1.  Significantly  differentially  expressed  genes 

1115  in  Flt3ITD/*  HSCs  and  HPCs.  Fold  changes  are  shown  for  270  probes  that 

1116  represent  254  differentially  expressed  genes  in  Flt3ITD/+  HPCs  relative  to  controls. 

1117  HSC  fold  change  data  are  also  shown.  P-values  and  adjusted  p-values  reflect 

1118  significance  across  the  entire  time  course. 

1119 
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1120  Figure  7  -  source  data  table  2.  Self-renewal-related  and  commitment- 

1121  related  gene  sets.  These  gene  sets  were  determined  by  identifying  genes  that 

1122  were  more  highly  expressed  in  HSCs  relative  to  HPCs  (self-renewal)  or  HPCs 

1123  relative  to  HSCs  (commitment).  In  each  case,  a  fold  change  threshold  of  5  and 

1124  an  adjusted  p-value  threshold  of  <0.05  were  used  to  define  the  lists. 

1125 

1126  Figure  8  -  source  data  table  1.  Heterochronic  gene  expression  in  wild  type 

1127  HSCs  and  HPCs.  Source  data  (log2)  are  shown  for  the  250  most  differentially 

1128  expressed  probes  in  developing  HSCs  (228  unique  genes). 

1129 

1130  Figure  9  -  source  data  table  1.  Flt3ITD  and  Runxl  mutations  cooperatively 

1131  induce  changes  in  gene  expression  in  post-natal  HPCs.  Genes  that  are 

1132  differentially  expressed  (Adjusted  p<0.05,  fold  change  >3)  in  Flt3ITD/+;  Runxl iV+ 

1133  HPCs  at  P0  and  PI 4  are  shown  in  the  indicated  tables. 

1134 
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Figure  5  -  figure  supplement  1 
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Figure  7  -  figure  supplement  1 
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Figure  8  -  figure  supplement  1 
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Figure  8  -  figure  supplement  2 


Euclidean  distances 
E14.5  to  P0<P0  to  Adult 
P14  to  Adult  <  E14.5  to  P14 
E14.5  to  PO  <  E14.5  to  P14 
P14  to  Adult  <  PO  to  Adult 


p  values 
<0.02 
<0.02 
<0.03 
<0.03 
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